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ABSTRACT 

PSRB1055— 52 is a middle-aged (r = 535 kyr) radio, X-ray, and 7-ray pulsar showing X-ray thermal 
emission from the neutron star (NS) surface. A candidate optical counterpart to PSRB1055— 52 was 
proposed by Mignani and coworkers based on Hubble Space Telescope (HST) observations performed 
in 1996, in one spectral band only. We report on HST observations of this field carried out in 2008, in 
four spectral bands. The astrometric and photometric analyses of these data confirm the identification 
of the proposed candidate as the pulsar's optical counterpart. Similarly to other middle-aged pulsars, 
its optical-UV spectrum can be described by the sum of a power-law (PLo) component (F u oc v^ a °) 1 
presumably emitted from the pulsar magnetosphere, and a Rayleigh- Jeans (RJ) component emitted 
from the NS surface. The spectral index of the PLq component, oq = 1.05 ± 0.34, is larger than 
for other pulsars with optical counterparts. The RJ component, with the brightness temperature 
To = (0.66 ± 0.10)d3 50 i?o 13 MK~ (where ^350 and Ro,i3 are the distance to the pulsar in units of 
350 pc and the radius of the emitting area in units of 13 km), shows a factor of 4 excess with respect 
to the extrapolation of the X-ray thermal component into the UV-optical. This hints that the RJ 
component is emitted from a larger, colder area, and suggests that the distance to the pulsar is smaller 
than previously thought. From the absolute astrometry of the HST images we measured the pulsar 
coordinates with a position accuracy of 0'.'15. From the comparison with previous observations we 
measured the pulsar proper motion, /i = 42 ± 5 mas yr _1 , which corresponds to a transverse velocity 
U t = (70±8)d 35 okms- 1 . 
Subject headings: pulsars: individual (PSRB1055— 52) — stars: neutron 



1. INTRODUCTION 

The radio pulsar B1055— 52 (also known as 
PSRJ1057-5226) was discovered by Vaughan & 
Large (1972). The period, P = 197 ms, and its 
derivative, P = 5.83 x 10~ 15 s s _1 , correspond to the 
spin-down age r = P/2P = 535 kyr, rotational energy 
loss rate (spin-down power) E — 3.0 x 10 34 erg s^ 1 , and 
surface magnetic field B = 1.1 x 10 12 G. Neither proper 
motion nor parallax measurements have been reported 
for this pulsar so far. Based on the pulsar's dispersion 
measure, DM=30.1 pc cm -3 , the distance was estimated 
as d = 1.53 kpc for the Galactic free electron density 
model by Taylor k Cordes (1993), while the NE2001 
model by Cordes & Lazio (2002) yields d = 0.73 ± 0.15 
kpc. 

PSRB1055— 52 was one of the few radio pulsars de- 
tected by the Einstein X-ray observatory (Cheng & 
Helfand 1983). Observations with EXOSAT showed 
a thermal nature of its soft X-ray emission, presum- 
ably originating from the neutron star (NS) surface 
(Brinkmann & Ogelman 1987), while ROSAT observa- 
tions found that at least two spectral components were 
needed to fit the spectrum (e.g., a blackbody [BB] and a 
power-law [PL]) and discovered X-ray pulsations at the 

radio period (Ogelman & Finley 1993). The similarities 
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in the ages and X-ray properties with two nearby middle- 
aged pulsars, PSRB0656+14 (r = 110 kyr) and Geminga 
(r = 340 kyr), for which ROSAT also detected thermal 
X-ray emission, earned these three pulsars the nickname 
of "The Three Musketeers" (Becker & Triimper 1997). 

The X-ray spectrum and pulsations of PSRB1055- 
52 were also studied with ASCA (Greiveldinger et al. 
1996; Wang et al. 1998) and BeppoSAX (Mineo et al. 
2002). Recent observation with the more sensitive Chan- 
dra (Pavlov et al. 2002) and XMM-Newton (De Luca 
et al. 2005) have shown that the phase-integrated X- 
ray spectrum, similar to those of the other two "Mus- 
keteers", can be fitted with three spectral components: 
a cold BB (T c ~ 0.8 MK), a hot BB (T H « 1.8 MK), 
and a PL component (photon index T w 1.7, luminosity 
L nonth _ 8 x io30rf2 50 er g s" 1 in the 0.5-10 keV band, 
where d^ 50 = d/(750pc)), presumably emitted from the 
pulsar magnetosphere. Phase-resolved spectroscopy of 
the Chandra and XMM-Newton data suggests that the 
changing projected emitting area of the hot component 
is responsible for the bulk of the phase variations (Pavlov 
et al. 2002; de Luca et al. 2005). 

PSRB1055— 52 was one of the seven 7-ray pulsars de- 
tected by the EGRET instrument aboard the Compton 
Gamma-ray Observatory (CGRO) (Fierro et al. 1993; 
Thompson et al. 1996). It has also been studied by 
the Fermi LAT detector (Abdo et al. 2010). Its 0.1-100 
GeV luminosity, I 7 ~2x lO 34 ^ 2 -^ erg s _1 (assuming an 
approximately isotropic emission), implies an extremely 
high 7-ray efficiency, n^ = L 1 jE ~ 0.6d^ 50 . 

To study and interpret the magnetospheric and ther- 
mal emission from PSRB1055— 52, the X-ray and 7-ray 
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data should be supplemented by optical data. How- 
ever, optical observations with ground-based telescopes 
(Manchester et al. 1978; Cheng & Hclfand 1983; Big- 
nami et al. 1988) failed to detect PSRB1055-52 because 
of the presence of a bright (V w 14) F-type field star 
w 4" from the radio pulsar position (dubbed Star A by 
Manchester et al. 1978). Thanks to the sharp resolution 
and high sensitivity of the Faint Object Camera (FOC) 
aboard the Hubble Space Telescope (HST), Mignani et 
al. (1997; hereafter M+97) were able to identify a faint 
(U ~ 25) probable counterpart at the radio pulsar po- 
sition. Since the observation was taken in only one fil- 
ter, the shape of the optical spectrum of the candidate 
counterpart remained unknown, and even its identifica- 
tion with the pulsar was not certain. To confirm the 
identification, measure the spectral shape (in particular, 
separate the magnetospheric and thermal components, 
similar to the two other Musketeers - see, e.g. Kargal- 
tev & Pavlov 2007; KP07 hereafter), and measure the 
optical pulsations, new HST observations were obviously 
required. Our proposal to reobserve PSRB1055— 52 with 
HST was accepted for Cycle 16 (program 11154). These 
observations are described in Section 2, while the results 
are presented in Section 3, discussed in Section 4, and 
summarized in Section 5. 

2. OBSERVATIONS AND DATA ANALYSIS 

As no timing capabilities were available in the HST Cy- 
cle 16, we could only image the target in several filters. 
To observe PSRB1055-52 in the far-UV (FUV), most 
useful for studying the thermal component, we employed 
the Solar Blind Channel (SBC) detector of the Advanced 
Camera for Surveys (ACS). The SBC is a Csl Multi- 
Anode MicroChannel Array (MAMA) photon-counting 
detector, which is operated in ACCUM mode producing 
time-integrated images. The detector is sensitive to ra- 
diation in the 1115-1700 A wavelength range, and it pro- 
vides a nominal field of view (FOV) of 34'.'6 x 30'.T, with a 
pixel scale of 0'.'034 x 0'.'030. To filter out the geocoronal 
lines, we used the longpass filter F140LP (pivot wave- 
length A = 1528 A; equivalent gaussian FWHM AA = 
297 A), which cuts off radiation at A < 1350 A. The 
pulsar was observed on 2008 February 13 over two con- 
secutive spacecraft orbits for a total integration time of 
5569 s. The data were reduced and calibrated (including 
the linearity correction, dark subtraction, flat-fielding) 
through the ACS data reduction pipeline (CALACS), using 
the closest-in-time calibration frames. Geometric distor- 
tion correction, cosmic-ray filtering, and exposure stack- 
ing were applied during the pipeline processing using the 
Multidrizzle task. This task also applied a re-sampling 
of the drizzled image to an even pixel size of 0'.'025. 

We also observed PSRB1055— 52 in three optical bands 
with the Wide Field and Planetary Camera 2 ( WFPC2) 4 
on 2008 March 8. To exploit the maximum spatial reso- 
lution, the pulsar was centred in the Planetary Camera 
(PC) chip, which has a pixel scale of 0'.'045 and a FOV 
of 35" x 35". To maximize the spectral coverage, we 
observed through the F450W (A = 4557 A; AA = 951 

4 We originally proposed the observations in the optical bands 
for the higher throughput ACS Wide Field Channel (WFC), but 
we had to switch to the WFPC2 since the WFC was turned off in 
January 2007 because of problems with the ACS electronics. 



A), F555W (A = 5442 A; AA = 1229 A), and F702W 
(A = 6917 A; AA = 1381 A) filters. Observations were 
carried out during four consecutive orbits for the total 
integration times of 3600 s in the F450W filter and 1800 
s in each of the F555W and F702W filters. For each 
filter, the observations were split in shorter, dithered ex- 
posures of 600 s to filter out cosmic ray hits and bad 
pixels. The data were processed through the WFPC2 
CALWP2 reduction pipeline for bias, dark, and flat-field 
correction and flux calibration. For each filter, we then 
combined single exposures with the IRAF task drizzle 
to produce co-added and cosmic-ray-free images. 

We also used the archival FOC observations by M+97. 
They were carried out on 1996 May 11 in the F342W 
filter (A = 3402 A; AA = 519 A) during three spacecraft 
orbits with the exposure times of 2771 s for the first orbit 
and 3062 s for the second and third orbits, corresponding 
to the total integration time of 8895 s. The data were 
obtained with the high-resolution F/96 relay, with the 
FOV of 7" x 7" and pixel scale of 07014. The data were 
retrieved from the HST science data archive 5 and on-the- 
fly recalibrated through the CALF0C reduction pipeline 
with the updated reference files. 

3. RESULTS 

3.1. Detection of the candidate counterpart 

The most recent published radio coordinates of 
PSRB1055-52 (Newton et al. 1981), 

a = 10 h 57 m 58?84 ± 0?03, 5 = -52°26'5673 ± 073, (1) 

are for the epoch of 1978.13 (equinox J2000). These co- 
ordinates, without proper motion values, are reported 
in the ATNF radio pulsar catalogue 6 (Manchester et al. 
2005), and they were also used by Weltevrede et al. 
(2010) to compute the PSRB1055— 52 radio ephemeris 
for the timing of the Fermi 7-ray observations. Since 
pulsars are known to have high velocities, PSRB1055— 52 
could move a few arcseconds in the 30 year span between 
the epochs of the radio position and our HST observa- 
tions (2008.12 and 2008.18). For instance, for an aver- 
age radio pulsar velocity of 400 km s _1 (e.g., Hobbs et 
al. 2005), we would expect a proper motion of 110 dj 50 
mas yr _1 and a position shift of 374 d? 50 in 30 years, in 
an unknown direction. Such an uncertainty hampers a 
direct identification based on the pulsar's position only. 
Fortunately, the presence of Star A in the immediate 
vicinity of the pulsar can be used for the pulsar identifi- 
cation. 

The 10" x 10" cutout of the SBC image of the 
PSRB1055-52 field (left panel of Figure 1) shows the 
only two objects detected in the entire SBC FOV, sepa- 
rated by a distance of 4741. From the comparison of this 
image with that obtained by M+97 with FOC, it is nat- 
ural to assume that the northwestern and southeastern 
objects are Star A and the candidate pulsar counterpart, 
respectively. 

Indeed, the nominal SBC coordinates of the northwest- 
ern source are a = 10 h 57 m 58?790 and 5 = -52°26'53"04. 
The most accurate coordinates of Star A, cxa — 

See http://archive.stsci.edu. 
6 See http://www.atnf.csiro.au/research/pulsar/psrcat/. 
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10 h 57 m 58 i :728±0!004 and 5 A = -52°26'52'.'45±0'/06, at 
the epoch of our HST observation, can be derived from 
the Third US Naval Observatory CCD Astrograph Cata- 
log (UCAC3; Zacharias et al. 2010), with account for the 
Star A's proper motion, \x a = —8.9 ±1.3 mas yr _1 and 
fig = 7.8 ± 2.2 mas yr -1 , given in the same catalog and 
independently verified by us (see Appendix). As the dif- 
ference between these and the nominal SBC coordinates 
(0'.'57 and —0759 in right ascension and declination, re- 
spectively) is within the error budget of the HST abso- 
lute astrometry, we conclude that the SBC northwestern 
source is indeed Star A. 

Applying the boresight correction to the SBC image, 
such that the coordinates of the northwestern source co- 
incide with the UCAC3 coordinates of Star A, we obtain 
the following coordinates of the southeastern source 



10 h 57 m 58!954, 



-52°26'56'/37, 



(2) 



with a nominal position uncertainty of about 0'.'07. 
This uncertainty is mostly determined by the uncer- 
tainty of the UCAC3 absolute astrometry, thanks to 
the small angular distance (4'.'41) between Star A and 
PSRB1055-52, and the - 0.1° accuracy on the HST 
roll angle with respect to the equatorial reference frame. 

The derived position is consistent with the candi- 
date pulsar counterpart position in the FOC observa- 
tion (a = 10 h 57 m 58?832, 5 = -52°26'56"28), within the 
uncertainty of the FOC position (w l'.'O) 7 and the un- 
certainty caused by an unknown proper motion since the 
epoch of the FOC observation (1996.36). This suggests 
that we indeed detected the PSRB1055— 52 candidate 
counterpart of M+97 in our SBC image, while the ex- 
tremely blue color (obvious from the comparison of the 
FUV (SBC) image with the optical (FOC and WFPC2) 
images; see Figure 1 and also Section 3.2) virtually proves 
that this object is indeed PSRB1055-52. 

Using the offsets between the pulsar counterpart and 
Star A measured in the SBC image (2702 and -3792 in 
right ascension and declination, respectively), we looked 
for the pulsar counterpart at the corresponding posi- 
tion in the WFPC2 images. We detected it in both the 
F702W image (seen in the circle labeled "SBC" in the 
right panel of Figure 1) and in the F555W image but not 
in the F450W one. 

3.2. Photometry 

We measured the magnitudes and fluxes of the 
PSRB1055-52 counterpart in the WFPC2, ACS/SBC, 
and FOC images through customized aperture photome- 
try using the IRAF package digiphot. For the WFPC2 
F555W and F702W bands, we used a small source aper- 
ture with the radius r = 3 pixels (07135) to maximize the 
signal-to- noise ratio (S/N) . We sampled the background 
in a surrounding annulus with the inner radius rt, = 5 
pixels and radial width 5r = 10 pixels; the distance of 
2 pixels from the source aperture was chosen to exclude 
the brighter portion of the point spread function (PSF) 

7 Before 2002, the pair of guide stars used to compute the as- 
tromctric solution in the HST focal plane were selected from the 
GSC1.0 (Lasker et al. 1990), which had a mean positional error 
of about l'.'O. After 2002, the guide stars arc selected from the 
GSC2.3 (Lasker et al. 2008), which has an improved mean posi- 
tional error of about 0'.'3. 



wings. We then applied the aperture correction using the 
fractional encircled energy (FEE) coefficients (Holtzman 
et al. 1995), and the correction to compensate for the 
time- and position-dependent charge transfer efficiency 
(CTE) losses of the WFPC2 detector (Dolphin 2009). 
For estimating the 3 a flux upper limit in the F450W 
band, we used the same r = 3 pixels source aperture 
but sampled the background in a wider annulus, Sr = 30 
pixels, around the expected pulsar position. 

For the SBC data, we chose a source aperture of r = 6 
pixels (072) and a background annulus with r\> = 20 pix- 
els and Sr = 10 pixels. The large radial distance of 14 
pixels between the source aperture and the background 
annulus was chosen to exclude the count excess over the 
PSF wings, which is seen from the comparison of the 
count distributions around the pulsar and around Star 
A (see the insets in the left panel of Figure 1). The 
3 a excess (e.g., 114 ± 38 counts in the 073-077 annulus) 
cannot be ascribed to different PSFs of the pulsar and 
Star A because the latter is only one magnitude brighter 
than the pulsar in the SBC F140LP band, and it is sep- 
arated from the pulsar by only 4741. From these data 
alone, we cannot exclude the possibility that the excess 
is a compact pulsar wind nebula (PWN), but, as the 
excess is not seen in other bands, it can also be just a 
small-scale enhancement of the detector background at 
the pulsar position. Since the brightness profile of Star 
A is closer to that of the model PSF, we used it as a ref- 
erence to compute the aperture correction by fitting its 
growth curve (i.e., the number of source counts within 
an aperture as a function of its radius). 

We also measured anew the magnitude and flux of the 
counterpart in the FOC image. Following the approach 
described by Pavlov et al. (1997), we chose an optimal 
source aperture from the measured growth curves for 
each of the three orbits. We found that the growth curves 
for the first and third orbits were consistent with each 
other and compatible with the nominal values of FEE 
at w 3400 A, while the growth curve for the second or- 
bit showed a peculiar behavior (e.g., no saturation up to 
r ~ 20 pixels, which means an unusually broad PSF). 
Therefore, we excluded the second orbit data from the 
analysis and chose r — 8 pixels (0711), r& = 10 pixels, 
and Sr = 10 pixels. 

In all cases, we applied the countrate-to-flux conver- 
sion by computing the photometric zero points in the 
STmag system using the image keywords PH0TFLAM and 
PH0TZPT, derived by the HST data reduction pipelines. 
To calculate the extinction-corrected fluxes, we esti- 
mated the interstellar reddening toward the pulsar, 
E(B — V) = 0.07, using the hydrogen column density 
N H = 2.7 x 10 20 cm" 2 , derived by De Luca et al. (2005) 
from the X-ray fits, and using the correlation between 
N H and E(B - V) found by Predehl & Schmitt (1995). 
We then calculated the extinction coefficients using the 
optical extinction curves by Fitzpatrick (1999) for the 
WFPC2 and FOC passbands and the UV extinction 
curves by Seaton (1979) for the SBC passband. The mag- 
nitudes and the measured, F° hs , and dereddened, F^ er , 
spectral fluxes are presented in the last two columns of 
Table 1. 

We plotted the dereddened spectral fluxes F^ er ver- 
sus frequency v in Figure 2. The shape of the broad- 
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band spectrum and, in particular, the relatively high 
FUV brightness of the object are inconsistent with it 
being a non-collapsed star in the Galaxy or an extra- 
galactic object. Therefore, the photometry proves un- 
equivocally that the source is the optical counterpart of 
PSRB1055-52. 

It is clear from Figure 2 that a simple one-component 
spectral model cannot fit the optical-UV spectral energy 
distribution (SED). However, similar to the optical-UV 
spectra of other middle-aged pulsars (e.g., KP07), the 
SED can be fitted by a model consisting of an (optical) 
power-law (PLo) component and a Rayleigh- Jeans (RJ) 
component: 



iC r = k 



VQ 



27r R l hT 2 



(3) 



where i?o is the radius of equivalent emitting sphere, 
and To is the brightness temperature (both as seen by 
a distant observer). At E(B - V) = 0.07, the fitting 
parameters are a = 1.05 ±0.34, f = 0.112 ±0.003 /zJy 
at i/ = 1 x 10 15 Hz, and R 2 T = (510 ± 81)d? 50 km 2 
MK (or T = (3.02 ± 0.48) df 50 Rq 2 13 MK > where #0,13 
is the radius in units of 13 km). The PLo component, 
which can be interpreted as magnetospheric radiation, 
dominates in the optical (A > 3000 A). The energy flux 
of this component is J-o° nth w 1.3 x 10 -15 erg cm -2 s _1 
in the 3000-9000 A band. The RJ component dominates 
in the FUV and is likely emitted from the NS surface. 
We will discuss implications of these results in Section 
4.1. 

3.3. Astrometry 
3.3.1. Absolute position of the pulsar 

The pulsar counterpart position given by Equation (2) 
was obtained from the SBC astromery using the Star A 
position in the UCAC3 catalog. Of course, using just 
one reference star may introduce a significant bias in the 
absolute target position. Fortunately, the WFPC2 has 
a much larger FOV than the SBC, allowing astrometry 
with larger numbers of reference stars. 

For the WFPC2 astrometry, we used two catalogs, 
GSC2.3 (Lasker et al. 2008) and 2MASS (Skrutskie et 
al. 2006). Although the astrometric accuracy of these 
catalogs is somewhat lower than that of UCAC3, they 
contain more stars per given area, which improves the 
statistics of the astrometric fits. Moreover, most UCAC3 
stars identified in the WFPC2 FOV are saturated, which 
prevents their use for the astrometric calibration. We 
used the F702W image, in which PSRB1055-52 was de- 
tected with a higher S/N than in the F555W image. We 
produced the WFPC2 mosaic image with the STSDAS 
task wmosaic, which also corrects for the geometric dis- 
tortion. Then, we measured the centroids of the reference 
stars (13 stars in GSC2.3 and 16 stars in 2MASS) in the 
WFPC2 pixel coordinates through gaussian fitting with 
the Graphical Astronomy and Image Analysis (GAIA) 
tool 8 and used the catalog sky coordinates of these stars 
to compute the pixel-to-sky coordinate transformation 
with the code ASTROM 9 . This yielded astrometric fits with 

8 See star-www.dur.ac.uk/~pdraper/gaia/gaia.html. 

9 See http://star-www.rl.ac.uk/Software/software.htm. 



the rms values a r — 0'.'17 for both the GSC2.3 and 
2MASS stars. To these we added in quadrature the un- 
certainties of the registration of the WFPC2 image on 
the catalog reference frames, Ctr = 0'.'14 and 0'/08, for 
the GSC2.3 and 2MASS, respectively. Following Lat- 
tanzi et al. (1997), these uncertainties were estimated 
as er t r = (3/ J V cat ) 1 / 2 CT cat , where 3 1 / 2 accounts for 3 free 
parameters in the astrometric fit, <7 ca t is the mean posi- 
tional error of the catalog coordinates (0'.'3 and 0'.'2 for 
the GSC2.3 and 2MASS, respectively), and A cat is the 
number of catalog stars used to compute the astrometric 
solution. Accounting for the 0'.'15 and 0'.'015 uncertain- 
ties on the link of the GSC2.3 and 2MASS to the Inter- 
national Celestial Reference Frame (ICRF), the overall 
uncertainties of the astrometry (1 a position errors) are 
S r = 0726 and 0719, respectively. 

We applied these astrometric solutions to compute the 
sky coordinates of PSRB1055— 52 from its pixel coordi- 
nates. As it has been done for the reference stars, we 
computed the pixel coordinates through a gaussian fit- 
ting, with an uncertainty of < 0701 for the pulsar cen- 
troid, negligible in comparison with the overall uncer- 
tainty of our absolute astrometry calibration. We then 
obtained the GSC2.3-based pulsar coordinates 

a = 10 h 57 m 58?961, 8 = -52°26'56717, (4) 

with the position error <5 r = 0726, and the 2MASS-based 
coordinates 



a = 10 h 57 m 58!967, 8 = -52°26'56731, 



(5) 



with 8 r = 0719. These two positions are consistent with 
each other and with the UCAC3-based position derived 
by using Star A as the only reference star in the SBC 
astrometry (Equation 2). 

Since the coordinates in Equations (4) and (5) have 
been computed using different catalogs and different ref- 
erence stars, they can be considered independent and 
hence can be averaged. The calculation of the weighted 
means yields 



10 h 57 m 58!965, <5 = -52°26'56"26, 



(6) 



with S r = 0715. Thus, the WFPC2 astrometry, re- 
calibrated with the GSC2.3 and 2MASS, has provided 
the absolute pulsar coordinates at the epoch of 2008.18, 
which is separated by about 30 years from the epoch of 
the only published radio position 10 . 

3.3.2. Proper motion of the pulsar 

We used the optical coordinates of PSRB1055— 52 to 
measure its proper motion from the comparison with 
the radio coordinates at the 1978.13 epoch (Equation 
1). We note that, since the GSC2.3 and 2MASS, as 
well as UCAC3, are linked to the ICRF, there is no a 
systematic offset between the optical and radio coordi- 
nates. Using Equation (6), we infer the displacement 
Aa cos 5 = +1714 ±0729 and AS = +0704±0732, which 
corresponds to the proper motion 

/j, a = +38 ± 10 masyr -1 , /j,g = +1 ± 11 masyr -1 (7) 

10 After completing our astrometric analysis, we became aware 
that an unpublished radio timing position was recently obtained 
(R. N. Manchester, private communication), which we found con- 
sistent with our own results. 
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in right ascenscion and declination. Using the UCAC3- 
based positions (Equation 2) results in virtually the same 
proper motion (/x Q = 35 ±9 mas yr" 1 , [j,$ = — 2± 10 mas 
yr -1 ), whose errors are dominated by the errors of the 
radio position (Equation 1). 

The pulsar proper motion can be independently mea- 
sured through relative HST astrometry, by comparing its 
position with respect to Star A in the SBC image (epoch 
2008.12) with that in the FOC image (epoch 1996.36). 
Of course, measuring a relative proper motion with only 
one reference star comes with caveats. Firstly, Star A 
is slightly saturated in the FOC image and part of its 
PSF is outside the FOC FOV (see Figure 1 of M+97), 
which lowers the precision of the star centroid determina- 
tion. To be very conservative, we assumed an uncertainty 
of 5 pixels (0'.'07). The pulsar centroids were measured 
with a much higher precision of < 0.1 pixels in both 
the SBC and FOC images (< 0"0025 and < O'.'OOW, re- 
spectively) , and the residual geometric distortion (about 
/ /004 and 0'.'007 rms for the SBC and FOC, respectively 
- see Apellaniz & Cox 2008 and Nota et al. 1996, respec- 
tively) is also much smaller than the FOC centroiding 
uncertainty for Star A. 

Secondly, with only one reference star, we have to rely 
upon the nominal values of the angles between the de- 
tector axes and the directions of right ascension and dec- 
lination, which may introduce an error in the direction 
of the proper motion. However, this error is very small 
in our case, thanks to the accuracy of the detector posi- 
tion angle and the small separation of Star A from the 
pulsar (see §3.1). Last but not least, the uncertainty of 
the proper motion of the reference star should be small 
enough not to hamper the measurement of the relative 
pulsar proper motion. We carefully verified it for Star A 
(see Appendix). 

Using the value of the pixel scale for the two detec- 
tors and the nominal orientations of the two images with 
respect to the equatorial reference frame, we found the 
displacement of the pulsar with respect to Star A during 
the time span of 11.76 years between the two HST obser- 
vations: Aa cos<5 = 0'/61±0'.'14 and AS = -0'.'19±0'/14. 
The uncertainty of the displacement is dominated by the 
large uncertainty of the Star A centroid in the FOC im- 
age. This displacement corresponds to the proper motion 
fi a = 52 ± 6 mas yr -1 and jj,s = —16 ± 6 mas yr . Cor- 
recting these values for the proper motion of star A as 
given in the UCAC3 catalog (see § 3.1), we obtained 



Ha = 43 ± 6 mas yr 



Us = —5 ± 6 mas yr 



(8) 



This proper motion agrees with that measured from the 
comparison of the radio and optical absolute coordinates. 
The weighted mean of the proper motion values deter- 
mined by the two methods is 



/j, a = 42 ± 5 mas yr 1 , /is = — 3 ± 5 mas yr 



-i 



or 



^ = 42 ±5 masyr -1 , P.A. = 94° ± 7 C 



(9) 



(10) 



for the total proper motion and position angle (counted 
east of north). The large proper motion provides addi- 
tional evidence that the M+07 candidate is indeed the 
pulsar counterpart. 



4. DISCUSSION 
4.1. Pulsar spectrum 

We have shown in Section 3.2 that the optical-UV SED 
for PSRB1055— 52 can be described by a sum of PLo 
(magnctosphcric) and RJ (thermal) spectra (see Equa- 
tion 3). It is interesting to compare these optical com- 
ponents with the corresponding spectral components at 
higher energies, particularly in soft X-rays where both 
the thermal and magnetospheric components are seen. 
As we have mentioned in Section 1, the pulsar's X-ray 
spectrum can be described by a model that consists of 
cold BB (BBc), hot BB (BB H ), and X-ray power-law 
(PLx) components. The BBc an d BBh components 
are presumably emitted from the bulk of the NS sur- 
face and polar caps, respectively, while the PLx com- 
ponent, F v = fx(E/Eo)~ a *, represents the magncto- 
sphcric X-ray emission. The parameters of the three 
components, as inferred from the XMM-Newton data by 
De Luca et al. (2005), are the following: T c = 0.79 ±0.03 
MK, R c = 12.3to 5 7 d 750 km; T H = 1.79 ± 0.06 MK, 
Ru = (0.46 ± 0.06) d 7 50 km; and a x = 0.7 ± 0.1, 

fx = L3i"o.i x 10 ~ 31 cr g cm ~ 2 s_1 Hz_1 at S = 1 
keV, respectively. These components and their extrapo- 
lations into the optical-UV domain are shown in the left 
panel of Figure 3. 

We see from Figure 3 that the extension of the best-fit 
PLx component into the optical overshoots the optical 
fluxes by a factor of ~ 3, similar to most of other pulsars 
detected in the optical (e.g., KP07; Mignani et al. 2010). 
However, in contrast to the majority of those pulsars, 
the PLo component is apparently steeper than the PLx 
component (ao = 1.05±0.34vs. ax = 0.7±0.1), and per- 
haps steeper than the PLo components in other pulsars 
(e.g., a = 0.41 ±0.08 and 0.46±0.12 in PSR B0656+14 
and Geminga, respectively; see KP07). Moreover, the ex- 
trapolation of the best-fit PLo component into the X-ray 
range docs not intersect the PLx spectrum, being at least 
a factor of ~ 10 lower at E > 1 keV. Such behavior is 
unusual (a notable exception is PSR B0540-69; Mignani 
et al. 2010), and it might suggest that different popula- 
tions of ultrarelativistic electrons are responsible for the 
optical and X-ray emission of PSRB1055— 52. However, 
taking into account the large uncertainties of the spectral 
slopes, it seems more plausible that ao is substantially 
smaller than its best-fit value (e.g., ao J$ 0.5), in which 
case the extrapolation of the PLo component can, in 
fact, be smoothly connected with the PLx spectrum. 

Interestingly, the slope a 7 = 0.06 ± 0.1 of the Fermi 
LAT spectrum, fit by a PL with an exponential cut- 
off (F v ex E- a - ■ cxp(- E/E cut ); Abdo ct al. 2010), is 
even flatter than the X-ray slope. On the other hand, 
as we see from the right panel of Figure 3, the optical 
and 7-ray points can be connected by a PL spectrum 
with the slope ao 7 ~ 0.46, which, however, goes slightly 
above the PL tail of the X-ray spectrum. We cannot rule 
out the possibility that further observations (especially 
in the IR-optical), supplemented by a joint multiwave- 
lcngth analysis, will show the magnetospheric spectrum 
of PSRB1055— 52 to be similar to those of other pulsars 
with optical, X-ray, and 7-ray counterparts, including 
the other two Musketeers. 

Another property of PSRB1055— 52 to compare with 
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those of other pulsars detected in the optical is the re- 
lationship between the optical and X-ray luminosities of 
the magnetospheric emission. For instance, for the wave- 
length range 4000-9000 A, the flux in the PLq compo- 



nent, Jq 



•nonth 



0.94 x 10 15 erg cm 2 s 1 , corresponds 



to the luminosity Lq 



nonth 



47rd 2 J"g onth = 6.3 x 10 28 d? 50 
erg s" 1 and optical efficiency r/o = L^^/E = 2.1 x 
10~ 6 <i 2 50 . At the adopted distance of 750 pc, both the 
luminosity and efficiency of PSRB1055— 52 are consider- 
ably higher than those of PSR B0656+14 and Geminga, 

4.2 x 10" 7 (i 2 



for which L^ onth = 1.6 x 10 28 d 2 



290 J ^O 



<290- 



and L^ onth = 4.8 x lO 27 ^, tj = 1.5 x lO" 7 ^, re- 
spectively (KP07). This hints that the actual distance 
to PSRB1055— 52 is smaller, which is supported by the 
analysis of the thermal emission below. On the other 
hand, the ratio of the optical (4000-9000 A) and X- 
ray (1-10 kcV) nonthermal luminosities, L™ nth /L™ nth = 
0.9 x 10~ 2 (which does not depend on the distance), is 
within the same range of ~ 10~ 3 -10 -2 as for all other 
pulsars with optical counterparts (Zavlin & Pavlov 2004; 
Zharikov et al. 2006). 

Let us now compare the optical-UV and X-ray ther- 
mal components of the PSRB1055— 52 spectrum. The 
BBh component is not seen in the optical because of 
the small emitting area, while the extrapolation of the 
BBc into the optical-UV is a factor of « 4 below the 
observed RJ component. Such an "optical-UV excess" 
of the RJ component has been seen in the nearby radio- 
quiet isolated NSs (RQINSs), such as RXJ1856. 5-3754 
and RX J0720.4-3125 (e.g., Kaplan 2008, and references 
therein), and in the millisecond pulsar J0437— 4715 (Kar- 
galtsev et al. 2004). It, however, has not been observed 
in middle-aged pulsars. In particular, the extrapolated 
BBc component lies slightly above the optical-UV RJ 
component in Geminga (Kargaltsev et al. 2005), while 
these components virtually coincide with each other in 
PSR B0656+14 (KP07). 

The nature of the optical-UV excess is not fully under- 
stood; a popular hypothesis is that it is due to a nonuni- 
form distribution of the temperature over the NS surface, 
such that the BBc component originates from a smaller 
and hotter region than the optical-UV RJ component 
(e.g., Pavlov et al. 2002). Such an assumption is crudely 
equivalent to adding a "very cold" thermal component, 
BBycj which dominates in the UV but makes very lit- 
tle contribution in X-rays. In other words, the product 
RqTo in the second term of Equation (3) is approxi- 
mately equal to RqTq + Ry C Tyc- Taking into account 

SlOc^o 



that RlT c « 119c&= n km 2 MK and R^Ti 



l c J c 



';-,{) 



km 2 MK, we obtain R? r ~T* 



vc-'ve 



391d? 50 km 2 MK (or 
Tyc ~ 2.31 d 2 50 i?y^ 13 MK). Unfortunately, we cannot 
measure Tyc independent of d /Ry C in the RJ regime. 
We can, however, estimate an upper limit on Tyc from 
the requirement that the BByc spectral flux extrapo- 
lated to soft X-ray energies becomes so small, in com- 
parison with the BBc, that it does not affect the X-ray 
fit. We found that it occurs at T vc < 0.45 MK (at this 
temperature the contribution of the BByc flux is < 10% 
above 0.3 keV - see Figure 3, left). 

Such a limit on Tyc implies a considerably smaller dis- 
tance to the pulsar than that estimated from the pulsar's 



dispersion measure. Indeed, from the above estimate on 
-Ryc^vc we obtain d < 25.4(i? vc /lkm)pc < 330i? N s,i3 
pc, where we took into account that i?vc is smaller than 
the NS radius, .Rns = 13-Rns,13 km. Since the NS ra- 
dius (as seen from infinity) can hardly exceed 20 km, the 
distance is < 500 pc, and it can even be substantially 
smaller than this upper limit. 

We should note that such estimates are, of course, 
model dependent. For instance, the temperature distri- 
bution over the NS surface can be smooth (i.e., there 
are no three distinct regions as we implicitly assumed 
above). Moreover, the spectrum of the thermal emission 
from a NS surface region can differ from the BB, which 
may lead to an excess of the actually emitted optical-UV 
thermal component over the extrapolation of the X-ray 
thermal component even if the surface temperature is 
uniform (e.g., in the case of a light-element, H or He, 
NS atmosphere; Pavlov et al. 1996). However, fitting the 
model light-clement atmosphere spectra to the observed 
X-ray spectrum always leads to a lower temperature and 
a larger value of R/d (i.e., to a smaller distance at a 
given NS radius) than those obtained from the BB fits 
(Pavlov et al. 1995). Therefore, it seems hard to avoid 
the conclusion on a smaller distance 11 , although its value 
cannot be determined accurately from the data available. 
We believe 200 pc and 500 pc are reasonable lower and 
upper limits, and we will scale the distance to 350 pc 
below. 

The conclusion that the distance is smaller than previ- 
ously thought has important implications. For instance, 
the radius of the BBc emitting region, Rq = 5.7^0.3^350 
km, becomes substantially smaller than a plausible NS 
radius. It means that T c = 0.79 ± 0.03 MK is not the 
temperature of the entire NS surface, and it should notbe 
used for comparisons with the NS cooling models. This 
removes the problem of PSR B1055— 52 being too hot and 
luminous (£c,boi = 4.2 x 10 32 d 2 50 erg s^ 1 ) for its age, in 
comparison with other middle-aged NSs and the predic- 
tions of standard NS cooling models (e.g., Yakovlev & 
Pethick 2004). It also means that PSRB1055-52 is not 
a "very slowly cooling NS" (as suggested by Kaminker et 
al. 2002), and there is no need to invoke an unusually low 
NS mass to explain its thermal emission (see Yakovlev & 
Pethick 2004 for references). To infer the true thermal 
luminosity and average surface temperature, and com- 
pare them with the predictions of the NS cooling mod- 
els, the parallax-based distance to the pulsar should be 
measured, and the distribution of the temperature over 
the NS surface should be determined from the phase- 
resolved spectral analysis of the pulsar's soft X-ray and 
FUV emission. 

Another implication of the smaller distance is that the 
magnetospheric luminosities and efficiencies in various 
spectral bands are lower than previously thought, being 
close to the luminosities and efficiencies of other pulsars. 
For instance, if d ~ 350 pc, its 7-ray efficiency, ?7 7 ~ 
0.13 d 2 35Q in the 0.1-100 GeV band (Abdo et al. 2010), 

11 Note a similar discrepancy between the DM-based distance 
estimate (750 pc, in the model by Taylor & Cordes 1993) and 
the distance determined from parallax measurements (288_ 27 pc; 
Brisken et al. 2003) for PSR B0656+14. Interestingly, Anderson 
et al. (1993) estimated the distance to be 280_ 50 pc using the NS 
atmosphere models by Shibanov et al. (1993) and ROSAT data. 
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although still rather high, becomes not so different from 
the efficiencies of other pulsars detected with Fermi LAT, 
and the optical efficiency, 770 = 4.6 x 10~ 7 <i3 50 , becomes 
comparable with those of PSRB0656+14 and Geminga. 

4.2. Pulsar astrometry and kinematics 

PSRB1055— 52 is the third rotation-powered pulsar, 
after Geminga (Caraveo et al. 1998) and PSRB0540-69 
(Mignani et al. 2010), for which an absolute position has 
been determined with a sub-arcsecond accuracy using 
optical astrometry techniques. Knowing the astrometric 
position will be useful for the precise timing analysis of 
the pulsar. 

The proper motion we measured for the pulsar's optical 
counterpart (Equation 11) corresponds to a tangential 
velocity V t — (70 ± 8)^350 km s . This value is well 
below the average 400 km s _1 for radio pulsars (Hobbs 
et al. 2005), but it is close to that of the Vela pulsar 
(V t ~ 65 km s" 1 ; Caraveo et al. 2001). 

The proper motion in Galactic coordinates, 



[ii = 39 ± 5 masyr , /if, = 15 ± 5 masyr 



(11) 



shows that the pulsar, whose current coordinates are 
I = 285.984° and b = +6.649°, is moving away from the 
Galactic plane. We attempted to use the proper motion 
to locate the pulsar's birth place, identify a parent star 
cluster (OB association), and estimate the actual age of 
the pulsar that can differ from its spin-down age. How- 
ever, to calculate the pulsar's trajectory back in time in 
the Galactic potential, one should know the present dis- 
tance (which is poorly known), the tangential velocity 
(proportional to the distance), and the radial velocity 
V r (which is unknown). Therefore, one has to calculate 
a large set of pulsar trajectories for various values of d 
and V r . Using the code developed by Vande Putte et 
al. (2010), we calculated the trajectories on a grid of 21 
distances (in the range of 200 to 1200 pc) and 21 radial 
velocities (—500 to +500 km s _1 ). We also selected a 
sample of young (< 100 Myr) clusters/associations (de 
Zeeuw et al. 1999; Dias et al. 2002) in a 25° radius cir- 
cle around the present position of the pulsar, extrapo- 
lated back in time their motion using the same code, 
and looked for closest approaches of the pulsar and clus- 
ter trajectories. The results, however, turned out to be 
rather ambiguous because different values of d and (es- 
pecially) V r resulted in quite different candidate parent 
associations and pulsar ages (e.g., we found 10 candi- 
dates for the age range of 400-700 kyr) , and the uncer- 
tainty of V t propagated back in time resulted in very 
large uncertainties on the minimum separations. These 
uncertainties becomes even larger if one also accounts for 
the errors on the cluster proper motions, distances, and 
radial velocities. To obtain a more certain solution, we 
have to wait for the measurement of the pulsar's parallax 
and more accurate measurements of the proper motion, 
together with observational constraints on the pulsar's 
radial velocity, which could come, e.g. from the mod- 
elling of a possible bow-shock PWN created in the ISM 
by the supersonically moving pulsar (see, e.g. Pellizza et 
al. 2005). 

5. SUMMARY 



Using the HST observations, we have confirmed that 
the candidate proposed by M+97 is indeed the optical- 
UV counterpart of PSRB1055— 52, the tenth rotation- 
powered pulsar with a secured optical identification (see 
Mignani 2010 for a recent review). From multi-band pho- 
tometry we found that, similar to the middle-aged pul- 
sars PSR B0656+14 and Geminga, its spectrum can be 
described by the combination of a PLo component, with 
the spectral index ao = 1.05+0.34, and a RJ component, 
with brightness temperature To = (0.66+0.10) d§ 50 i?Q 2 13 
MK. The PLo component is steeper than those of 
PSRB0656+14 and Geminga and, possibly, of all the 
other rotation-powered pulsars with known optical coun- 
terparts. Moreover, it is possibly steeper than the X-ray 
PL component, which might suggest that different pop- 
ulations of relativistic electrons are responsible for the 
X-ray and optical-UV magnetospheric emission. To ver- 
ify this assumption, IR observations of PSRB1055— 52 
would be particularly useful, supplemented by a joint 
spectral fit of the IR-optical, X-ray, and 7-ray compo- 
nents, using multiwavelength spectral models for magne- 
tospheric pulsar emission. 

The observed RJ component exceeds by a factor of «4 
the extrapolation of the X-ray thermal component into 
the UV-optical, unlike PSRB0656+14 and Geminga but 
similar to RQINSs. It indicates that the temperature 
distribution over the NS surface is nonuniform, and the 
X-ray thermal component comes from a relatively small 
hotter region, while the main contribution to the optical- 
UV RJ component comes from a larger and colder area. 
It implies that the distance to the pulsar is considerably 
smaller than estimated from the pulsar's dispersion mea- 
sure and the models for Galactic electron distribution. It 
also means that the conclusions inferred from the com- 
parison of the X-ray temperature and thermal luminosity 
of PSRB1055— 52 with the NS cooling models require a 
revision. New HST observations in the UV and near- 
IR, including phase-resolved spectroscopy, are required 
to separate the thermal and magnetospheric components 
and infer the temperature distribution, while the mea- 
surement of the radio parallax is crucial to determine 
the actual distance and the NS radius. 

By recalibrating the astrometry of the HST images 
with the GSC2.3, 2MASS and UCAC3 catalogs, we 
have measured the PSRB1055— 52 absolute coordinates 
with a radial position accuracy of 0'/15 at the epoch of 
2008.18, about 30 years after the measurement of the 
only published pulsar's coordinates. From relative as- 
trometry between the HST images taken 12 years apart, 
we obtained the first measurement of the pulsar's proper 
motion, /1 = 42 ± 5 mas yr _1 with a position angle of 
94° ± 7°, which corresponds to the transverse velocity 
V t = (70 ± 8)d 350 km s" 1 . Further HST and radio ob- 
servations will allow one to measure the proper motion 
with a higher accuracy, to determine the distance from 
the parallax measurement, and find the true pulsar age 
and the parent cluster/association. 

Support for Program number 11154 was provided by 
NASA through a grant from the Space Telescope Science 
Institute, which is operated by the Association of Uni- 
versities for Research in Astronomy, Incorporated, under 
NASA contract NAS5-26555. RPM thanks A. C. Jackson 



Mignani, Pavlov, & Kargaltsev 



for her contribution in the reduction of the VLT data. 
We thank R. N. Manchester for communicating the radio 
timing position of PSRB1055— 52 before publication, A. 
Dolphin for useful discussions of the corrections for CTE 
losses in WFPC2, D. Vande Putte for making available 
to us his code for Galactic orbit simulations, and M. van 
Kerkwijk for his valuable advice on the search for the 



pulsar's birthplace. We also thank the anonymous ref- 
eree for the very careful reading of the manuscript and 
the useful suggestions that helped us to improve the pre- 
sentation of the results. The work by GGP and OK was 
partially supported by NASA grant NNX09AC84G. 
Facilities: HST (ACS, WFPC2),VLT (FORS1) 



APPENDIX 
STAR A PROPER MOTION 

Since the magnitude of the Star A proper motion is crucial for determining the pulsar's absolute position and proper 
motion (see Sections 3.1 and 3.3.1), we verified the UCAC3 value (/i Q = —8.9 ± 1.3 mas yr _1 ; [is = 7.8 ± 2.2 mas 
yr _1 ) against other astrometric catalogs. Unfortunately, star A is too faint (V « 14.6) to be in the Hipparcos catalog 
(V < 12.4; Perryman et al. 1997), while it is strangely not in the deeper Tycho-2 (V < 15.2; H0g et al. 2000) and 
UCAC2 (R < 16; Zacharias et al. 2004). It is, however, in the USNO-B1.0 catalog (B < 22; Monet et al. 2003), which 
reports a proper motion of fi a = +30 ± 4 mas yr _1 and fig = ± 6 mas yr _1 (much larger than in UCAC3), and 
in the PPMXL catalog (Roeser et al. 2010), which gives fi a = —13.8 ± 8.2 mas yr _1 and us = +4.5 ± 8.2 mas yr _1 , 
consistent with the UCAC3 value (albeit less accurate). 

As an independent test, we measured the relative proper motion of Star A. We used archival Very Large Telescope 
(VLT) images of the PSRB1055-52 field taken with FOcal Reducer/low dispersion Spectrograph ( FORS1) on 2002 
March 21, which yields a time baseline of about 6 years with respect to the WFPC2 observations. The images were 
taken in high resolution mode (O'.'l/pixel) but with a windowing of the CCD (1.'6 x lf6 FOV). A sequence of eighteen 
140 s dithered exposures was obtained through the Bcsscl B filter (A = 4290A; A A = 880 A). We applied standard CCD 
reduction with the MIDAS package and performed the astrometric calibration as described in Section 3.3 (position 
error S r = 0'.'3). Unfortunately, the image quality of the FORS1 data (seeing ~ l'/2) does not allow one to detect the 
pulsar in the halo of Star A. 

We then compared the relative position of Star A measured in the FORS1 (2002.21) and WFPC2 (2008.18) images. 
We used the F450W WFPC2 image, where Star A is not saturated and the pivot wavelength is similar to that of 
the Bessel B filter. We registered the WFPC2 image onto the FORS1 image, aligned along right ascension and 
declination to better than ~ 0.05°, by computing the transformation between the pixel coordinates of 15 stars seen 
in both images. This yielded an rms of 0'.'015 along each of the two axes, larger than the object centroid uncertainty 
(0'.'003). We thus determined the offset of Star A to be Aa cos S = —2.2 ± 15 mas and AS = 24.6 ± 15 mas, which 
corresponds to the 3cr upper limit on the proper motion of ~ 15 mas yr _1 , consistent with both the UCAC3 and 
PPMXL values. 
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TABLE 1 

HST PHOTOMETRY OF THE PSRB1055 — 52 COUNTERPART. 



Filter 


log V 


t 


r 


C r " 


m r 


Am r 


AniCTE 


wiob B 


pobs 


pdor 




[Hz] 


s 


pix 


counts 










A*Jy 


A*Jy 


WFPC2 F450W 


14.818 (0.045) 


600 


3 


<13.7 


>25.61 


-0.27 


-0.37 


>24.97 


< 0.26 


< 0.34 


WFPC2 F555W 


14.741 (0.049) 


600 


3 


27.3 ±5.0 


25.90 ±0.20 


-0.24 


-0.23 


25.43 


0.241 ± 0.044 


0.295 


WFPC2 F702W 


14.637 (0.043) 


600 


3 


30.3 ±3.5 


26.46 ±0.13 


-0.26 


-0.13 


26.08 


0.214 ± 0.024 


0.249 


SBC F140LP 


15.293 (0.042) 


5569 


6 


525 ± 24 


22.83 ± 0.05 


-0.19 




22.64 


0.248 ±0.011 


0.417 


FOC F342W 


14.945 (0.033) 


2916 


8 


242 ± 13 


25.18 ±0.06 


-0.18 




25.00 


0.140 ± 0.008 


0.190 



Note. — Columns report the camera and the filter, the pivot frequency and the band width (in parentheses) in logarithmic units (as derived 
from the SYNPHOT tables), the integration time of the average-combined frames, the photometry aperture radius r in detector pixels, the total 
number of counts (C r ) measured within the aperture r, the corresponding magnitude in the STmag system (m r ), the aperture correction (Am r ), 
the CTE correction (AmcTE), the observed magnitude m b s , and the observed (F° fa ) and extinction-corrected (F^ CI ) spectral fluxes. 



a For the WFPC2 and FOC observations, counts have been measured on the average-combined frames, 
average of the first and second exposure (see § 3.2). 



For the FOC, we used only the 








Fig. 1.— SBC F140LP (FUV band; left panel) and WFPC2 F702W (red band; right panel) images of the PSRB1055-52 field (10" X 10"). 
North is to the top, East to the left. The two objects in the SBC image are the pulsar (in the smaller thin circle of 0"3 radius, corresponding 
to a nominal uncertainty of SBC position) and Star A (4'.'4 northwest of the pulsar). Zoomed images of the pulsar and Star A, demonstrating 
the count distributions in and around these sources, are shown in the insets; the radii of the circles in the insets are 0'.'3 and 0"7. The 
small thin circle in the WFPC2 panel, labeled "SBC" is at the same position with respect to Star A as in the left (SBC) panel. In each of 
the panels the smaller thick circle of 0"45 radius is the error circle around the radio pulsar position at the epoch of 1978.13 (Newton et al. 
1981), while the larger thick circle of 3" 4 radius encloses possible positions of the pulsar at the epoch of our HST observations if the pulsar 
transverse velocity (in an unknown direction) were 400 km s — 1 , at d = 750 pc. The thin circle labelled "FOC" corresponds to the position 
of the PSRB1055— 52 candidate counterpart of M+97; its radius of l'.'O corresponds to the uncertainty of a nominal FOC position. See 
§3.1 for more detail. Notice the much higher relative brightness of the pulsar in the SBC image. 



HST observations of PSRB1055-52 
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Fig. 2. — Optical-UV spectrum of PSRB1055— 52. The spectral fluxes have been corrected for interstellar extinction assuming E(B — V) = 
0.07. The dashed lines correspond to the Rayleigh- Jeans (RJ; red) and power-law (PLo; blue) components of the best fit, the solid line 
shows the sum of these components. The shaded areas show lcr uncertainties of the fit. 
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Fig. 3. — Left: Comparison of the optical-UV and X-ray spectral models for PSRB1055— 52. The X-ray fit for the BBc + BB H + 
PLx model is taken from De Luca et al. (2005). The excess of the RJ component over the extrapolation of the X-ray BBc component 
into the optical-UV range could be explained by the presence of a "very cold" thermal component, BBvo The dotted blue curve shows 
this component for the maximum allowed temperature Tvc = 0.45 MK. The shaded areas show la incertainties of the fits. Right: 
Multiwavelength SED for PSRB1055-52, from optical to 7-rays, with spectral models for separate ranges. The 7-ray points are from 
Guillemot (2009), the spectral model from Abdo et al. (2010). The line v~ °- 46 shows a PL spectrum approximately connecting the optical 
and 7-ray bands. See Section 4.1 for more details. 



